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Ginzburg Landau theory for d-wave pairing 
and fourfold symmetric vortex core structure 
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The Ginzburg Landau theory for _y2-wa,ve superconductors is constructed, by starting from 
the Gor'kov equation with including correction terms up to the next order of ln(rc/r). Some 
of the non-local correction terms are found to break the cylindrical symmetry and lead to the 
fourfold symmetric core structure, reflecting the internal degree of freedom in the pair potential. 
Using this extended Ginzburg Landau theory, we investigate the fourfold symmetric structure 
of the pair potential, current and magnetic field around an isolated single vortex, and clarify 
^ l' concretely how the vortex core structure deviates from the cylindrical symmetry in the d.j.2_y2- 

D ' wave superconductors. 

c/^ ■ 
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J> §1. Introduction 

, Much attention has been focused on a vortex structure in high- Tc superconductors. It is expected to be clarified 
' how the vortex structure of the high- Tc superconductors is different from that of the conventional superconductors. 
0^ : By a number of experimental and theoretical investigations, it is concluded that the symmetry of these supercon- 
' ductors is most likely to be da;2_j,2-wave. Therefore, one of the points is how the vortex structure of d2^2_y2-wave 
. superconductors is different from that of isotropic s-wave superconductors. The d^2 _y2-wd.'ve pairing, i.e., — ky in 
I momentum space, has fourfold symmetry for the rotation about the c-axis. We expect that, reflecting this symmetry, 
the core structure of an isolated single vortex may break the cylindrical symmetry and show fourfold symmetry in 
(ia;2_j,2-wave superconductors. p. 
^ . In their experiments, Keimer et alu reported an oblique lattice by a small-angle neutron scattering study of the 
O vortex lattice in YBa2Cu307 in a magnetic field region of 0.5T< H <5T applied parallel to the c-axis. The vortex 
^ : lattice has an angle of 73° between the two primitive vectors and is oriented such that the nearest-neighbor direction 
Y\ ' of vortices makes an angle of 45° with th£|_a-axis. The oblique lattice was also observed by scanning tunneling 
J> I microscopy (STM) by Maggio- Aprile tt alu u They also observed the elliptic-shaped STM image of the vortex core, 
and concluded that this oblique lattice cannot be explained by considering only the effect of the intrinsic in-plane 
anisotropy, that is, the difference of the coherence lengths between a-axis and &-axis directions.&CF It is suggested 
' that this deformation from the triangular lattice in d^2_j^2-wave superconductors is due to the effect of the fourfold 
symmetric vortex core structure.BB'B 

The fourfold symmetric vortex core structure in (i^2_j^2-wave superconductors was so ^ix decived theoretically when 
the s-wave component is induced around a vortex of (i^2 —^i-wave superconductivityfl'oEHj'lla' This mixing scenario 
was mainly studied by Bcrlinsky alui and Ren et aZUeP in the two-component Ginzburg Landau (GL) theory 
for s- and d-wave superconductivity. According to the consideration based on the two-component GL theory, it is 
possible that the s-wave component is coupled with the d-wave component through the gradient terms. Therefore, 
the s-wave component may be induced when the-id-wave order parameter spatially varies, such as near the vortex 
or interface under certain restricted conditions The induced s-wave component around the vortex is fourfold 
symmetric. The resulting vortex structure in cZ-wave superconductors, therefore, exhibits fourfold symmetry. 

For this scenario to be applied to high-Tc superconductors, the amplitude of the induced s-wave component 
should be comparable to the already-existing main c?3.2_j^2-wave component. The_amplitude of the induced s-wave 
component is proportional to the s-wave component Vs of the pairing interactionJiH' Then, in the case where \ys\ is 
negligibly small compared with the dominant c?^2_j,2-wave pairing interaction, the induced s-wave order parameter 
is negligible. In this paper, we consider this pure dr^2_y2-^SNQ case, that is, the case when the induced s-wave 
component can be neglected. In our opinion, as the first step to understanding the fourfold symmetric vortex core 
in da;2_j,2-wave superconductors, the pure dj.2_j,2-wave case should be studied before considering the mixing of the 
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induced s-wave component. In this case, we have to consider a new scenario to explain the fourfold symmetric vortex 
core structure. 

When the induced s-wave component can be neglected, the order parameter reduces to be one component, i.e., only 
the da;2_j,2-wave component, and the GL equation reduces to the same form as that of the isotropic s-wave case within 
the conventional GL theory. Within this framework, the vortex core structure remains cylindrically symmetric, and 
the_vortex lattice forms a triangular lattice, even in the d3,2_j^2-wave superconductors. On the other hand, Ichioka et 
alx3 showed that the vortex core structure is fourfold symmetric in the pure d^2_y2-wave superconductors by using 
the quasi-classical Eilenberger theory. This means that we have to modify the conventional GL theory to explain 
the fourfold symmetric vortex core structure. For this purpose, we consider the following scenario. 

Strictly speaking, the conventional GL equation is valid only near the transition temperature Tc- Far from Tc, 
we have to include several correction terms of the order ln(Tc/T), which consist of both the higher powers of the 
order parameter and higher order derivatives, i.e., the so-called nonlocal terms. Among them, some nonlocal terms 
break cylindrical symmetry in the dx^-y2-wa.ve pairing case. Therefore, these correction terms lead to the fourfold 
symmetry of the vortex structure. 

So far, these terms neglected in the conventional GL theory were considered to explain the fourfold symmetric 
behavior of the upper critical field Hc2 in the a6-planell3' and the deformation of the vortex lattice from a triangular 
lattice. As for the fourfold symmetric behavior of HfS^ Takanaka and Kuboyao^ showed it from the GL theory 
with nonlocal correction terms, and Won and MakiQ'H from the Gor'kov equation. A stable vortex lattice was 
discussed by Won and Maki in the pure c?^2_j,2-wave superconductors. As for the case of the isotropic s-wave 
pairing but anisotropic density of states at the Fermi surface, Takanaka and Nagashimallj^ derived the GL equation 
with including the higher order correction terms, and studied the anisotropy of Hc2 and a stable vortex lattice. 
However, the anisotropic vortex core structure has not been studied so far. In this paper, we discuss the fourfold 
symmetric vortex core structure by considering the effect of these correction terms. 

The purposes of this paper are, first, to construct the GL theory from the Gor'kov equation for the pure d^^^y^- 
wave superconductors with including correction terms of the order ln(Tc/T) and, second, to investigate the fourfold 
symmetric core structure (i.e., pair potential, current and magnetic field) of an isolated single vortex by using this 
extended GL theory. We clarify how the vortex core structure deviates from the cylindrical symmetry. Here we 
consider the case of an isolated single vortex under a magnetic field applied parallel to the c-axis (or z-axis) in the 
clean limit. The Fermi surface is assumed to be two-dimensional, which is appropriate to high- Tc superconductors, 
and isotropic in order to clarify effects of the d-wave nature of the pair potential on the vortex core structure. The 
additional anisotropy coming from, e.g., the Fermi surface can be incorporated into our extended GL framework. 

The rest of this paper is organized as follows. In we construct the GL theory with including correction terms 
for the dx2-y2--wsLve pairing. By using this GL theory, we investigate the fourfold symmetric vortex core structure. 
The pair potential is studied in and current and magnetic field in The summary and discussions are given in 
§|[ We set ?i = c = fee = 1 throughout the paper. 

§2. GL theory for d^2_y2-waLve pairing 
We consider the GL theory in pure dx2-y2-weLve superconductors. The pair potential is given by 

A{r,k):^T]{r)V2coa2e. (2.1) 

Here, r = (x, y) = (r cos (p, r sin (j)) is the center of mass coordinate and fe the relative coordinate of the Cooper 
pairs. Now, k is denoted by an angle 9 measured from the a-axis (or x-axis) in the a5-plane. To study the effect of 
correction terms in the order of small parameter ln(Tc/r) ~ 1 — T/Tc, we consider the GL theory by including the 
next higher order terms of rj and its derivatives other than those in the conventional GL theory. The GL equation 
and the current density j are derived as follows from the Gor'kov equation within the weak coupling approximation 
(see the detailed derivation in Appendix A), 



In ( ^ ) ,7 - /? ( ryhp - ^£)^,; 
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(higher order terms) — 0, (2-2) 



jx = 2joIm 
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where 
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= 0.908. 
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21C(3) 
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49C(3)2 

with Riemann's ^-functions C(3) and C(5), the Fermi velocity vp and the density of states at the Fermi surface 
Nf- The differential operator D (D*) is defined by £) = V + i{2n/(j)o)A {D* = V - i(27r/(^o)^) with the vector 



potential A and the flux quantum (pQ. In Eqs. (2.2) and (2_^), D'^D''^ means to take all the permutations of the 
product, such as D^Dy = D^Dy 



D^DyD, 



It is convenient to discuss Eqs. (2.2) and (2.3) in the following dimcnsionless unit 



DyDl 



V, 



with the GL coherent length ^ = {f3vj,/6\n{Tc/T)y/'^ and the energy gap 770 = {\n{Tc/T)/py/'^. 
In the dimcnsionless unit, the GL equation in Eq.(2.2) and j in Eq.(O), respectively, are written as 



(2.5) 



ri-r]\ri\'^+D'^rj + alrL I ^ 



-i{4h|2(r>2^) + TjHD'i^r + 2rj\D^\^ + SiDTj^v*} 



+ 0({ln(T,/T)}2) = 0, 



(2.6) 



= 2Im 



+ {D^Tj){Dlf^r + {D^Dlr^)rf + {Dyfj){D,Dyrjr] + 0({ln(T,/T)}2) 



(2.7a) 
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l-3|rypaln 



{Dyr^W + ^In ('^] {7{Dlrj)rf - 7{Dl^){Dyr,r 



-{Dyrj){Dlr^r + {DlDyi^W + {D,ii){D,Dyr,r] + 0({ln(r,/r)}2 



(2.7b) 



In Eqs. (2.6) and (2/7), the neglected terms, which are higher order of rj and derivative, are in the order {ln(rc/r)}" 
(n > 2). We note that, in the limit T Tc, Eqs. (2.6) and (2.7) reduce to those of the conventional GL theory, 
and are the same forms as those of the conventional s-wave case. This indicates that there are no difference between 
s-wave and d-wave pairing within the conventional GL framework. Therefore, the vortex structure for d-wave pairing 
remains cylindrically symmetric within the conventional GL framework. The difference between s-wave and d-wave 
pairing first appears in the correction terms of the order ln(Tc/r). In particular, the ZJ^-terms in the GL equation 
(2.6) and the Z?'^-terms in the current density (2.7) are seen to break the cylindrical symmetry, and lead to the 
fourfold symmetry. These correction terms play an important role as the temperature decreases below the transition 
temperature. 

§3. Pair potential around a single vortex 

3.1 Symmetry consideration 

The (ia;2_j,2-wave pairing, i.e., k^ — ky = cos 20 in momentum space, has fourfold symmetry for the rotation about 
the c-axis. Reflecting this property, the symmetry of a single vortex in d^-2_y2-^&Ye superconductors is in the class 
D^^\D2) X R. In the limit T Tc, the GL equation (2.6) reduces to the conventional GL equation, which is the 
same ap-.that of the isotropic s-wave pairing case. Therefore, following the well known consideration about a single 
vortex,EJ) we obtain the structure of a single vortex as follows in the limit T ^ Tc, 



A(r,fc) = r){r)V2cos2e = /o(r)e*"*%/2 cos 26*, 



(3.1) 
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which is the isotropic state with the winding number 1 around the vortex. As noted by Volovik,0 eight elements 



of the symmetry operations of the function A(r, fe) in Eq. (3.1) form the group D4{E) 



Di{E) =. {E, C2e'\ U2^^R, U2^ye'^R, C^e'^^^ Cle'"^'', V^.^+ye'^ '^R, U2,^-ye~''' '^R}- (3.2) 

Here, R is the time-reversal, C„ the (27r/n)-rotation about z-axis, and Vi.xi C^2,t/7 t^2,x+y and U2,x-y denote the 
TT-rotation about the x-axis, the y-axis, the lines x + y = and x — y = 0, respectively. 

On lowering temperature, other components with different winding numbers may be induced in ri{r). If the 
symmetry D^{E) is conserved, the possible pair potential is generally restricted to the following form (brief derivation 
is given in Appendix B), 

oo 

'7(0 = (3-3) 

n— — OQ 

where /n(r) is real. In the following, we determine the amplitude /n(r) from the GL equation (2.6). 

3. 2 Pair potential around a single vortex 

We determine the amplitude /n(r) in Eq.(3.3) up to the order ln(Tc/r). Therefore we expand /„(r) in the powers 
of ln(Tc/T) as follows, 

Ur) = fl^\r) +aHTjT)fi'Hr) + 0({ln(T,/r)}2), (3.4) 

where fn\r) — for n ^ since A(r, k) reduces to Eq. ( |3.l| ) in the limit T Tc- 

We notice that, as we consider the isolated single-vortex in the extreme type II superconductors (GL parameter 
K ^ 1), the vector potential in D can be neglected. Ej-* 



Substituting Eq. (3.4) into Eq. (2.6) and writing the differential operators by the cylindrical coordinate, we obtain 



the following simultaneous differential equations, 

(r) ^ Cir) + V^(l)/(")(r) - C (rf = 0, (3.5) 

F!,'\r) ^f^'\r) + V\l)f^'\r) - 3/^(0^^'^^ + J^i(l)/o^°^(0 

-I (5/r (r)v^(i)/r (.) + 5/r (.)/r'('f - ^.i'\rf) + ^/r = o, (3.6) 

Fi'\r) ^ f['\r)+V^5)f['\r) - (2/{i)(r) + /iV(r)) (r)^ + (!)/(") (r) = 0, (3.7) 

F!V(r) ^ /i\)(r) + V2(-3)/i\)(r) - (2/i\\r) + f['\r)) C {r)' + 1^0^(1) f^^\r) ^ 0, (3.8) 

F(i)(r) ^ /(i)(r) + V\4n + l)fi'\r) - (2/(i)(r) + /«(r)) f^'^r)' = 0, (n ^ 0, ±1) (3.9) 
where the differential operators V^(n), Di{n) and D±{n) are defined as 

VN = TT + -1 2' (3-10) 

dr^ r ar 

^ , , 2 d^ I d^ Id , / 2 d2 2d 4 \ 71^ 

- ^ + - + - -^^, + ^)+^' (3-11) 

^ _d^ 6 d3 ^ 15 d^ 15 d / I d^ 6 d^ 14 d 12 

dr^ r dr-^ dr^ r*^ dr \ r dr^ 
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As for the boundary condition, 



/r(0) = /l^'(0) = (3.13) 
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at the vortex center so that divergence does not occur at r 



0, and 



oo 



5 

12' 



3) = (n 7^ 0) 



(3.14) 



far from the vortex core, where the r-dependence of ^^(t) can be neglected ensuring rjir) to become the bulk value. 
Equation ([3.14 ) is obtained by setting the terms with the differential operators to be in Eqs. ( |3.5| ) - (^). 
Equation (3^) corresponds to the conventional GL equation, leading to the well known vortex core structure of 



the conventional GL theory for Jq {r First, we obtain /g (r) from Eq. (^.5[), and substitute it to Eqs. (3.6) 



Then, we obtain the correction terms fQ^\r) from Eq. (3.6), and f±i{r) from Eqs. (3/7) and (li). From Eq. 
(3.9), we obtain fn\r) = for n 7^ 0,±1. In the order ln(Tc/T), therefore, the components with e~'^"^ and e^*"^ are 
induced in addition to the isotropic components of e*"^. Components with other winding numbers are induced in the 
higher order of ln(Tc/r), for example, components with e~^"^ and e^"^ first appear in the order {ln(Tc/T)}^. 

^''^<r) in Eqs. 



The concrete form of the solution /. 
we set 



q) - (3.9) cannot obtained by the r- or r ^-expansion. When 



m—1 



to study the solution for r ^ 1, we obtain the fohowing results for m < 3, 



(0) 
0,m 



C^ol = 0, C%^ = 0, (m^3) C« =0 for 0,-1, 



(3.15) 



(3.16) 



which means that /o(r) = 0(r), f-i{r) = 0{r^) and /i(r) = O(r^). However, the other coefficients „j cannot be 
determined uniquely in this r-expansion. On the other hand, when we set 



(3.17) 



to study the solution for r 3> 1, inconsistency among Eqs. (3.5) - (3.9) occurs in determining c[^^_^ and Ci\^ It 



suggests that the singular terms such as r ' e ^ probably exist in addition to r '"-terms. 
3.3 Numerical Results 



We solve numerically the differential equations (3.5) - (3.£) by using the so-called relaxation method, which is 
often used in determining the vortex structure from the GL theory. The relaxation step 



jO")(^)[now] ^ /^j)(r)[oid] _^ cF^^\r) (c : constant) 



(3.18) 



is repeated until the condition \F^ ''\r)\ ^ 1 is fulfilled for each r, where F.^'" (r) is defined in Eqs. (| 

Figure 1 (a) shows fQ^\r), which is the correction to the factor of the isotropic component e*"*. As seen from Eq. 
(3.4), the correction ahi{Tc/T) f^^^r) is added to the solution of the conventional GL theory flf\r). In Fig. 1 (b). 



(3.9) 



f^{r) and /i\''(r) are presented, where /i\''(r) behaves as near the core while f[^\r) as r^. These components 
induce the fourfold symmetric structure of the pair potential. They both attain a maximum at several coherence 
length from the origin, and smoothly and slowly fall off. It is to be noted that they never drop exponentially. It 
indicates that the effect of the fourfold symmetric structure becomes small but remains far from the vortex core, 
and may affect a vortex-vortex interaction and the formation of a stable vortex lattice. This point may become 
important when considering a stable vortex lattice and its orientation relative to the underlying crystal in d-wave 
superconductors. It is confirmed that the other components fn^\r) {\n\ > 2) reduce to in the relaxation method. 

Using these solutions, the pair potential around a single vortex is written as 

V{r) = {Cir) +aHT,/T) (/(^'(r) + fi'\r)e^^'^ + /iV(r)e-4'*) } + 0({ln(T,/r)}2). (3.19) 
Then, the amplitude of the pair potential around a single vortex is given by 

\v{r)\ = [{C ir) + « HTc/T) {/^^^ (r) + cos 40 (/f ^ (r) + (r)) } }' 
+ {a ln(T,/r) sin 40 (/f ^ (r) - /i^ (r)) }'] 
= (r) + aln{T,./T) [f^'\r) + cos 40 (/f ^(r) + /«(r)) } + 0({ln(r,/T)}2). (3.20) 
As for the amplitude |77(r)| in Eq. (3.20), we see that the deviation from the cylindrical structure behaves like cos 40 
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Fig. 1. The spatial variation of the components of the order parameter around a vortex in the dimensionless unit: (a) /g^^(r'), (b) 
/-['^'(r) (solid curve) and /i\'('r) (long dashed curve), associated with the phase factor e"^, e^"^ and e~^"^, respectively. 
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Fig. 2. The spatial profile of the order parameter around a vortex, (a) The deviation of the amplitude |»7(r, (/))| at = 0° from |r?(r, 0)| 
at </) = 45° along the radial direction, where the factor of aln(rc/T) in Eq. (3.21) is displayed, (b) The correction of the phase, 
*(r, 4>), at = 22.5°, where the factor of aln(Tc/T) in Eq. (3.23) is displayed. 



around a vortex in the order \vJ(rclT). In Fig. 2 (a), we show the r-dependence of the difference of the amphtude 
along the 0° direction (along the x-axis and y-axis) and along the 45° direction (along the line y — ±a;), i.e., 



|ry(r,0 = 0)| - W,^ = 7r/4)| = 2aln(Tjr) (/f >(r) + /^(r)) . 



(3.21) 



In the figure, the factor of Q;ln(rc/r), i.e., 2(/|^-'(r) + /i\''(r) ) is plotted. It indicates that the anisotropy localizes 
around the core region, and that the amplitude along the 0° direction is suppressed compared with that along the 
45° direction. 



As for the phase of the pair potential, we consider the correction of the phase to the leading term which goes like 
gi0_ Then, we introduce <I>(r) defined by 



ry(r) = |77(r)|e"^+'*('^). 



(3.22) 



From Eq.(3.19), the correction of the phase is given by 



GL theory for d-wave pairing and fourfold symmetric vortex core structure 



7 



^[r) = tan 



MTc/T)sm^[fi'\r)-f^^l{r)) 



Cir) +aHn/T) {f^'\r) + cos 40 (^f['\r) + /iV(r)) }^ 
aln(r,/r) sin 40 [fi'\r) - /«(r)) //f (r) + 0({ln(r,/r)}2). 



(3.23) 



It indicates that $(r) behaves hke sin 40 around a vortex in the order \n(Tc/T). In Fig. 2 (b), we show the 
r-dependence of <I>(r) along the 22.5° direction from the a axis, i.e., $(r, = tt/8). In the figure, the factor of 
aln(Tc/T), i.e., (/f'(r) - /^l {r) ) / f^°\r) is plotted. It is seen that the anisotropy of the phase rather extends 
compared to that of the amplitude in Fig. 2 (a) , and that $(r) > for the region < < 7r/4. 

The above-mentioned behavior of the pair potential r]{r) qualitatively agrees well with that obtained in the 
quasi-classical Eilenberger theory (see Figs. 1 and 2 in Ref. 12 for comparison). 

§4. Current and magnetic field around a vortex 

The current density and the magne tic fi eld aroun d a single vortex is considered with the correction terms up to 
the order hi{Tc/T). Substituting Eq. (3.19) to Eq. (2.7), we obtained the radial component, jr — jx cos -I- jy sin 0, 
and the rotational component, = —j^ sin + jy cos 0, of the current density as follows. 



>(r) - aln(re/T)sin40(ji(r) + j2{r)) + OUHn/T)}^), 
Mr) = ^/rW +«Mr./T){-j3(r) + cos 40(^-1 (r) + j4(r)) } + 0({ln(T,/T)}2), 



(4.1) 
(4.2) 



where 



j,{r) ^ 2 (/r(^)/f ^'(^) - /r(r-)/i\^'(r) + /iV(r)/r'(r) - f['\r)C\r)) , 



(4.3) 
(4.4) 



Mr) - ^/r - ^/r'M/r w + 1 (/r'w^ - 2/r"w/r m + - j/r 'm/^^^'w, (4.5) 



U{r)^^.i"\r) (3/«(r)-/«(r)) 



(4.6) 



\j{r)\ = Idr)^ + alniTjT)[-Mr) + cos40(ji(r) + Mr)) } + 0({ln(T,/T)}2), (4.7) 



From Eqs. (4.1) and (4.2), the amplitude of the current density is given by 
2 

r " 

whic h is the same expression as in Eq. ( [4.2|) within the order ln{Tc/T). In the limit T Tc, Eqs. ( |4.l| ), (|4. 
and (17) give the cylindrically symmetric structure of the conventional GL theory. This profile of |j(r)| is shown 
in Fig. 3 (a) as a function of r. 

From Eq. (4/7), we see that the deviation of |j(''')l from the cylindrical structure first appears in the order ln(rc/T) 
and behaves like cos 40 around a vortex. In Fig. 3 (b), we show the r-dependence of the difference of \j{r)\ along 
the 0° direction and along the 45° direction, i.e.. 



iir, = 0)1 - b-(r,0 = 7r/A)\ = 2a ln(r,/T) (ji (r) + j4(r)) . 



(4.8) 



In the figure, the factor of a\n{Tc/T), i.e., 2(ji(r) + j4{r) ) is plotted. It is seen that \j{r, = 0)| < \j{r, = 7r/4)| 
near the vortex core and \j{r, = 0)| > \j{r, = 7r/4)| far from the core. It should be noted that this sign change 
in Fig. 3 (b) occurs at a point where /i\ (r) in Fig. 1 (b) reaches the minimum. 

The resulting magnetic field 6(r) can be found from the Maxwell equation: 

j{r) = K^V X 6(r) (4.9) 

in the dimensionless unit: Eq. ( |2.5D and 

b{r)/V2HcK^b{r), (4.10) 
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Fig. 3. The spatial profile of the current density around a vortex, (a) The amplitude |j(r)| as a function of r in the limit T —> 
(b) The deviation of \j{r,(l>)\ at </) = 0° from \ j{x,(j>)\ at </> = 45° as a function of r, where the factor of aln(Tc/T) in Eq. (4.i 
displayed. This plot shows the deviation of |j(r)| from the cylindrical structure. It is noted that the sign of the deviation changes at 
r = 3.2 from negative to positive. 



where He — 1/k^^2-\/2 is the thermodynamic critical field. From Eqs. ( [4.lD , ( |4.2| ) and (4.9), the magnetic field, 
which has only the z-component, is given by 

h,{r) =b,{r = 0)-^!^ lflfi\r' fdr' + -l^aln{TjT)[{l - cos4(/.)r(ji(r) + j2(r)) 

-^l^\jiir')-j3{r')+Mr'))dr']+0{{HTjT)y). (4.11) 

In the limit T Tc, Eq. (4.11) gives the cylindrically symmetric structure of the conventional GL theory. This 
profile of bz{r) — bz{r = 0) is shown in Fig. 4 (a) as a function of r. 

From Eq. (4.11), we find that the deviation of bz{r) from the cylindrical structure first appears in the order 
ln(Tc/r) and behaves like cos 40 around a vortex. In Fig. 4 (b), we show the r-dependence of the difference of bz{r) 
along the 0° direction and along the 45° direction, i.e., 

bz{r,(b = 0) - fe.(r,0 - V4) = -^aln(re/r)r(ji(r) +j2(r)). (4.12) 

In the figure, the factor of K~^a\n{Tc/T), i.e., —r{ji{r) + j2{'i') )/2 is plotted. It indicates that bz{r) is large along 
the 0° direction compared with along the 45° direction, that is, the magnetic field extends along the 0° direction 
around a vortex. It is noted that, refiecting the sign change in Fig. 3 (b), the anisotropy of bz{r) becomes weaker 
far from the vortex core. 



The above-mentioned behaviors of j{r) and bz{r) qualitatively agree well with those obtained in the quasi-classical 
Eilenberger theory (see Figs. 4 - 7 in Ref. |l2|for comparison). 

§5. Summary and discussions 

First, the GL theory is constructed from the Gor'kov equation with including the correction terms up to the order 
ln(rc/T) in pure dx2_y2--wave superconductors. In the limit T Tc, the GL equation and the current density in 
the (ia;2„j,2-wave case reduce to the same form as that of the isotropic s-wave superconductors. Thus, within the 
conventional GL framework, the vortex core structure remains cylindrically symmetric even in the (i^2„j,2-wave case. 
The difference between the s-wave and d^2_y2--wsLve first appears in the correction terms of the order ln(rc/r). 

Second, by using this extended GL theory, we investigate the fourfold symmetric vortex core structure in the 
da;2_j,2-wave superconductors. Among the correction terms, there are some terms which break the cylindrical sym- 
metry and have fourfold symmetry, which is the fact reflecting the fourfold symmetry of the Cooper pair, — ky. 
Then, the vortex core structure becomes fourfold symmetric. We study the fourfold symmetric structure of the pair 
potential, the current and the magnetic field around an isolated single vortex with the correction terms up to the 
order ln(Tc/T). As a result, the amplitude of the pair potential is suppressed along the 0° direction compared with 
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Fig. 4. The spatial profile of the magnetic field around a vortex, (a) The magnetic field (r) — bz{r = 0) as a function of r in the 
limit T Tc, where the factor of in Eq. (4.f 1) is di splay ed, (b) The deviation of bz{r, <f>) <j> = 0° from bz{r, (f>) at (f> = 45° as a 
function of r, where the factor of K~^oln(Tc/T) in Eq. (1.12) is displayed. This plot shows the deviation of bz{r) from the cylindrical 
structure. Far from the vortex core, this deviation decreases and the fourfold symmetric structure becomes weaker. 



along the 45° direction. The phase of the pair potential also shows fourfold symmetry around a vortex. Along the 
45° direction, the amplitude of the current is enhanced near the vortex core and suppressed far away from the vortex 
core compared with along the 0° direction. Reflecting this behavior of the current, the magnetic field is enhanced 
along the 0° direction near the vortex core, and this fourfold symmetry of the magnetic field fades out slowly far 
from the vortex core. These fourfold symmetric structures of the pair potential, the current and the magnetic field 
qualitatively agree well with the results by the quasi-classical Eilenberger theory.H On lowering temperature, the 
correction terms of the order ln(rc/r) make important roles, and the fourfold symmetric vortex core structure 
becomes clear. 

It should be noted that, while our obtained vortex core structure has similar fourfold symmetry to that predicted 
by Xu et alia from the two-component GL equation, the origin of the fourfold symmetry is quite different from 
theirs. In their theory, the fourfold symmetry is induced by the mixing of the s-wave component. Then, in the case 
induced s-wave component is negligibly small, i.e., the pure d^2_y2--wsLve case, the vortex core structure reduces to 
cylindrically symmetric one. On the other hand, we consider the pure d^2_y2--w8ive case. Even in this case, the 
vortex core structure exhibits fourfold symmetry in our theory. It is due to the fourfold symmetric terms which are 
higher order of \n{Tc/T) and neglected in the conventional GL thepry. p. j— . 

In connection with these correction terms, Takanaka and KuboyatJ) and Won and MakiU'EJ' explained the fourfold 
symmetric behavior of Hc2 in the a6-pane and the deformation of the vortex lattice from a triangular lattice by 
considering the effect of these terms. 

By the extended GL theory with the correction terms, we succeed in reproducing qualitatively the fourfold 
symmetric vortex core structure obtained by the quasi-classical Eilenberger theory.ll3 This fact indicates that the 
correction terms of the order ln(r/rc) are important in the GL framework, when we consider the vortex structure 
in pure dx2-y2-wave superconductors. The GL framework presented here is found to be useful to analyze the vortex 
core structure and to obtain further physical insight. While enormous numerical computations are needed in the 
Eilenberger theory, the GL theory gives the information about the vortex structure with concise calculations. 

As shown in Fig. 1 (b), the effect of the fourfold symmetric vortex core structure becomes small but remains far 
from the vortex core. Then it may affect a vortex-vortex interaction and the formation of a stable vortex lattice. This 
point may become important when considering a stable vortex lattice and its orientation relative to the underlying 
crystal in d-wave superconductors, which belongs to our future study. 
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Appendix A: Derivation of the GL equation 



The GL equation (2.2) and the current expression (2.3) in the dx2_y2-wsLve pairing case are derived from the 



Gor'kov equation with including both the higher power of the pair potential and higher order derivative. Here, 
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we consider the singlet anisotropic pairing case, in general. Similar derivation was presented by Takanaka_and 
Nagashima in the case of isotropic s-wave pairing but anisotropic density of states at the Fermi surface.tS In 
our derivation, higher order terms are considered up to the order D™'Tf (m + n < 5), and the integral ^ dr and 
/ fifc/(27r)^ are noted by and Y^j^, respectively. 

As for the anisotropic pairing, the pair potential has the non-local form A(r, r'). Thus, in the Gor'kov equation for 
the Green function G^^^ (r, r') and th£.anomalous Green functions F^^^ (r, r') and F^^ (r, r'), the term A{r)F^^^ (r, r') 
in the isotropic s-wave pairing casecS is replaced by X^r" ^(''"i ''"'O^wn (''""i ^'^d the self-consistent condition is 
modified to 

A* (r, r') = V{r - r')T ^ F^^ (r, r'). (A-l) 

From the Gor'kov equation, we obtain the relations 

G.„(r,r') = GO Jr,r')- GZjr,r,)A{r,,r[)Fl(r[,r'), (A-2) 

ri,r[ 



Fl{r,r') = J2 G°^Jri,r)A*(ri,r;)G.„K,r'), (A-3) 

ri,r[ 

by using the Green function of the normal state, 

G°„(r-,r') = 5:-^e'^'-(^-^'). (A-4) 



For convenience, we introduce the center of mass coordinate f = (r + r') /2 and the relative coordinate r = r — r' 
of the Cooper pair. Then, we have r — f + f/2, r' = r — r/2 and Vr — Vr' = 2V:. The pair potential and the 
pairing interaction are, respectively, assumed to be 

A(f, k) = Y^ e-'^-^A(r, r') = ii{r)c^{k), (A-5) 



Vik -k') = Y^ e-'('^-'^')-V(f ) = V(i,*{k)(j){k'), 



(A-6) 



where (t){k) is a symmetry function, for example, <f){k) — '\/2cos2^ in the (i^2_j,2-wave pairing case and 4'{k) — 1 in 
the isotropic s-wave pairing case. Then, the Fourier transformation of A(r,r') is given by 



A(r,r') = ^e*''-+*'^-^(g)</'(fc). 

q,k 



From Eqs. (A-5) and (A-6), the self-consistent condition ( |A-1| ) is rewritten as 

i^*[r)^VTY,Y.^{k)Y, e'^-'Fl (r , r ') 

t^n k f 

On the other hand, the current density is given by 

j{r) - ^(V. - V.OT^G.„(r,r') 



2i\e\ 



T5^V = G.„(r,r') 



(A-7) 



(A-, 



(A-9) 



r=0 



From Eqs. ( A-2 ) and ( [A-3| ), Guj„ (t, r') in Eq. ( A-9 ) and F^ (r, r') in Eq. ([A-^ ) are expanded in the power of A as 
follows. 



where 



(r-, r') = G0„ (r, r') + G^^^' (r, r') + G^ (r, r') + - • • , 
i^i = F}}''^ (r, r') -f (r, r') -f F}}^^ (r, r') + - - • , 



=- E G[i„(r,ri)G0^„(r2,rl)G^Jr^,r')A(ri,r;)A*(r2,r'2), 

E (^2, rl)GO„ (r^, r3)G°_„ (r4, r!,)G:i„ {r'„ r') 



(A-10) 
(A-11) 

(A-12) 
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X A(ri, r;)A*(r2, r^)A(r3, r^)A*(r4, r^), (A-13) 

F^i'Hr,r') =- G^„(ri,r)G°„(r;,r2)G^Jr3,r',)G°„(r^„r') 

xA*(ri,r;)A(r2,r'2)A*(r3,rJ,), (A-15) 

xA*{r,,r[)A{r2,r'^)A*{r,,r',)A{r,,r'^)A*ir,,r',). (A-IG) 



We substitute Eqs. (|A-4|) and (|A-7] ) into Eqs. (A-12) - (A-16), and perform the integrals in Eqs. (|A-8|) - (A-16). 
Since q corresponds to the spatial variation of the pair potential (|q| ^ 1/^) and fc the relative motion of the Cooper 
pair (|fc| kp), that is, |q| ^ |A;|, we can take 

£k+q/2^£k + -z-k-q, (l){k + q) ~ (j){k), (A-17) 



where fc = fc/|fc|, and Eqs. (A-8) - (A-16) are calculated by expanded also in the power of k ■ q. The fc-integral is 
treated as J^k ~ I / d9/2TT. We write J d6/2Ti{- ■ ■) = {■ ■ for simplicity. 
As a result, the GL equation ([A-8|) and the current density ([A-9|) are written as 



ry*(r) =YVNpYA^{^{-ir\m?).^,nq)e~^1- 

q n—Q 

oo 

+ Y^Np J2 ({2^n+2 + «l+^3Mn+3}(-^ir(-^2)"H-^3)"M-^4n0(fc)r/^ 

Qi'^Qz ni~n4— 

(n— ni+n2+n3+n4) 

+ Y 6^^F^4(|0(fc)r)-r,*(qi)r;(q2),7*(q3)r7(q4)r;*(q5)e^(-^i+'?^-'?3+'/4-'/5)-r^ (A-IS) 



^1,92 rii, 712=0 

+ J2 mNFVFA,(^k{2i, + i^ + i, + 2i,)\m\')fJl{qiW{q2)v{q^W (A-19) 
where = wpfc • and ^„ is defined as follows: 



In^, (forn = 0) 

^"="^ ^ Tfc^(l-2^)^(" + l)' (for even n> 2) (^-20) 



Ci^Ti = — oo 



0, (for odd n) 

with the Debye frequency LUp, and the Euler constant 7. In the second term of Eq. (A-18), we define q^ = q^ —92 + 93- 
In the last terms of Eqs. (A-18) and (A-19), we present only the contributions of the lowest order of q. 

By performing the q-integral in Eqs. (A-18) and (A-19), q^viQi)^^'^^'^ are replaced by (— iD,i)"?7(r,i)|r;=r , where 
Di = Vr; + i{2Tr/(l)o)A{ri) is introduced instead of Vr^ for conserving the Gauge invariance. Then (A-18) and 
(A-19) are rewritten as 



77*(r) ^NpVj2An{HP*r\m\^)j^V*{r) 



00 

+ NpV Y ({2A„+2+*(K+^;Mn+3}(-*|'n"H^I>2)"^(-*|>;)"n^|)4)"1<^(fc)|^^^ 

rii ~n4— 
(n—ni +712 +713+714) 
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X77*(ri)?7(r2)77*(r3) 



+ 6NFVAM<l>{k)f)\7j{r)\'r,*{r), 



(A-21) 



+4z|e|7V;^i;i.A4(fc(-2|)i+^;-^3 + 2^:)|0(fc)|4^^r,(ri)?7*(r2)??(r3)ry*(r4) 



ri=r2=r3=r4=r 



(A-22) 



where ^ = vpk ■ D. In Eq. (A-21), we define D4 = -Dl - D2 - D^. 

We substitute 0(fe) = \/2cos20 (in the dx2_y2-wa.ve case) and k ■ D ^ DxCosd + Dy sin 0, and perform the 
^-integral. Then, using the relation {NpV)'^ = ln{2uj d7 / t^Tc) , we ob tain Eqs. (2.2) and (pT^). 

We confirm that the GL equation (2.2) and the current expression (2.3) are also derived by expanding the quaai- 
classical Eilenberger equation in the powers of D and 77 following the method presented by Schopohl and Tewordt.Ej-* 



Appendix B: Derivation of Eq. ( |3.3| ) 

Around a vortex, in general, r]{r) of the pair potential A(r, k) in Eq. (Ol) can be expanded by e*'"* (/: integer) 



as follows, 



Then, as for the U2,xR and U2,x-ye "^^i? in Eq. ( [3.2[ ), the transformations of A(r, k) are given by 

C/2,.ii'A(r,fc) = ^/;(r)e^'^V2cos20, 
I 

U2,x-ye-'''^^RA{r,k) = - ^ /;*(r)e-^('+i)^/2+*'"* ^2 cos 26*. 



(B-l) 



(B-2) 



(B-3) 



If the symmetry D4^{E) is conserved, Eqs. ( |B-2| ) and (B^) should be equal to A(r, k). Then, we obtain the relation. 



.fiir) = nir), 



-i(; + l)7r/2 _ 



= -1. 



(B-4) 



From Eq. (B-4), //(r) is real, and I = An + 1 {n: integer). Writing /;=4n+i(r) as /n(r), we obtain Eq. (3.3) 
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